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ABSTRACT 

Recent particle-in-cell simulations suggest that a large fraction of the energy dissipated in a 
relativistic shock is deposited into a Maxwellian distribution of electrons that is connected 
to the high-energy power-law tail. Here, we explore the observational implications of such 
a mixed thermal-nonthermal particle distribution for the afterglow and prompt emission of 
gamma-ray bursts. When the Maxwellian component dominates the energy budget, the after- 
glow lightcurves show a very steep decline phase followed by a more shallow decay when the 
characteristic synchrotron frequency crosses the observed band. The steep decay appears in 
the X-rays at ~100 sec after the burst and is accompanied by a characteristic hard-soft-hard 
spectral evolution that has been observed in a large number of early afterglows. If internal 
shocks produce a similar mixed electron distribution, a bump is expected at the synchrotron 
peak of the v fy spectrum. 
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1 INTRODUCTION 

After the prompt emission phase is over, the relativistic ejecta re- 
sponsible for the gamma-ray burst (GRB) drive a relativistic shock 
into the circumburst medium. The afterglow emission that follows 
the GRB is believed to come from the shocked external medium. 
Electrons are thought to be accelerated at the shock front and ra- 
diate throughout the electromagnetic spectrum via synchrotron and 
synchrotron-self-Compton (SSC) mechanisms (see Piran 2005 for 
a review). 

Particle acceleration in relativistic shocks is poorly under- 
stood. Most of the theoretical studies of GRB afterglows param- 
eterize the distribution of particles downstream of the shock by 
assuming that a fraction of the energy dissipated at the shock 
goes to accelerating all electrons into a pure power-law distribu- 
tion over a wide range of energies (Paczynski & Rhoads 1993; 
Sari et al. 1998). The possibility that only a fraction of the 
electrons are accelerated and the rest thermalize with temperature 
= kT/m^c^ ~ r, where T is the bulk Lorentz factor of the shocked 
plasma, has also been considered (Eichler & Waxman 2005). 

Recent particle-in-cell (PIC) simulations have begun to shed 
light on the acceleration processes at work in relativistic shocks. 
The picture that emerges for the particle distribution downstream 
from the shock is somewhat different than what has been tradition- 
ally assumed. Diffusive shock acceleration does appear to operate 
in such shocks, accelerating a small fraction of particles (few %) 
into a power-law distribution that carries some ~ 10% of the dis- 
sipated energy. The rest of the electrons are, however, found to be 
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heated to temperature ~ Tnip/mg carrying most of the dissipated 
energy (Spitkovsky 2008a,b; Martins et al. 2009). 

In this work, we explore the observational implications of such 
a mixed thermal-nonthermal particle distribution. We show how the 
synchrotron emission spectrum and lightcurves are modified com- 
pared to conventional afterglow models. We suggest that the spec- 
tral and temporal properties of the early X-ray lightcurves that show 
very steep decay may be produced with such a particle distribution. 
We also briefly discuss the possible implications of a mixed elec- 
tron distribution on the prompt GRB emission. 



2 MODELING THE ELECTRON DISTRIBUTION 
DOWNSTREAM FROM THE SHOCK 

The GRB afterglow emission is believed to come from the shock 
driven by the ejecta into the circumburst medium. Since the after- 
glow is steadily observed over a wide range of frequencies ranging 
from the radio to X-rays (Costa et al. 1997; Frail et al. 1997; van 
Paradijs et al. 1997; and likely extends into the y-rays), a broad 
particle distribution downstream of the shock is needed to explain 
observations (Sari et al. 1998; see also Paczynski & Rhoads 1993; 
Waxman 1997; Wijers et al. 1997). A power-law distribution in en- 
ergy for the particle number oc y^'' with index p ~ 2 - 2.5 is 
typically inferred (e.g., Waxman 1997; Galama et al. 1998). Par- 
ticle acceleration through repeated shock crossings is a promising 
mechanism to accelerate particles to such a distribution (e.g., Gal- 
lant, Achterberg & Kirk 1999). 

Consider a relativistic shock which has a relative Lorentz fac- 
tor F » 1 between the shocked and unshocked fluid. A rather eco- 
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Figure 1. The electron-number distribution for different values of the frac- 
tion S of the total energy that goes into the power-law component. For 5=1, 
the pure power-law distribution is reproduced (dashed line). The most re- 
cent PIC simulations indicate that 5 ~ 0. 1 (solid, black line). For 5 < 1/3, 
the Maxwellian component in the electron distribution is pronounced. 



nomical, and popular, parameterization of the particle distribution 
downstream of the shock assumes that all the electrons are accel- 
erated into a power-law distribution for y > y„„„. When p > 2, 
most of the energy of the electrons resides at the low end of the 
distribution. Then, 7,„,„ is found from the total available energy per 
particle assuming that a fraction of the energy dissipated by the 
shock goes into the electrons, resulting in y,,,,,, ~ e^Tnip/mg. In this 
parameterization, the distribution has a sharp cutoff below 7m,„. 

Eichler & Waxman (2005) considered the possibility of only a 
fraction of the electrons being accelerated to a nonthermal distribu- 
tion. The rest were assumed to isotropize and thermalize at a char- 
acteristic temperature ~ F given by the upstream kinetic energy 
of the electrons with no extra heating from the ions. In the Eich- 
ler & Waxman (2005) approach, the particle distribution is peaked 
at 7 ~ r and y ~ Tnip/nii, with most of the energy carried by the 
nonthermal component. The nonthermal electrons in this model are 
responsible for most of the emission (with the possible exception 
of the early radio afterglow). 

Recent PIC simulations of relativistic shocks have revealed 
both the stochastically accelerated and Maxwellian components in 
the downstream electron distribution (see Spitkovsky [2008b] and 
Sironi & Spitkovsky [2009] for simulations of pair shocks, and 
Spitkovsky [2008a] and Martins et al. [2009] for electron-ion sim- 
ulations). Simulations of relativistic shocks in weakly magnetized 
electron-ion plasma show that the thermalized electrons are sub- 
stantially heated with respect to their upstream bulk flow energy, 
receiving a large fraction of the dissipated energy from the ions. 
The electron thermal distribution is smoothly connected to the low 
end of the power-law component (see Fig. 1). The thermal com- 
ponent is a robust prediction of all simulations, and represents the 
downstream state of the bulk flow. It can be depleted if the trans- 
fer of energy to the power-law component is very efficient. Current 
simulations show that in shocks that accelerate particles, the power 
law tail can take of the order of 10 - 20% of the flow energy. As the 
present simulations necessarily run for a limited amount of time, it 
is possible that the steady state has not been achieved yet. There- 
fore, it makes sense to explore the observational signatures of the 
mixed thermal-nonthermal distribution for a range of acceleration 
efficiencies. 



In this work, we define the normalized electron distribution 
downstream of the shock by a continuous function connecting a 
relativistic Maxwellian to a power-law at y = y,,,;, : 



C<'(7,0), for 7<7,„;„ 



(1) 



CNf(y„„,&)(y/y„,„)- 



for 



r > 7 ml; 



where N'J'iy,®) = y^ exp(-y/&)/2&^ is the Maxwell distribution 
in the » 1 limit (which is of interest here) and C is a normal- 
ization constant. The ratio 7,„;,/0 controls the fraction S of the total 
energy that resides in the nonthermal component: 



£jNA7,®)dy 
f^yNAy,e)dy' 



(2) 



Furthermore, we assume that a fraction of the energy dissi- 
pated in the shock is picked up by the electrons. The average (ran- 
dom) Lorentz factor per particle is 



< 7 >E 



yNAy, &)dy ■ 



(3) 



Using the last expression with eqs. l[T), l|2j the electron distribution 
is determined (i.e., and 7,,,;, can be calculated) after choosing 6, 
€g, F and p. 

The mixed thermal-nonthermal distribution reduces to the 
pure power-law for 6 = 1. In this case 7,„,„ = 7,,,;, = (p - 
2)eeTmpl{m,{p - I)) which defines the minimum cutolf of the dis- 
tribution. In the limit of negligible nonthermal component (5 <k 1, 

=< 7 > /3 = e^Fmp/Bm,,. 



3 SYNCHROTRON SPECTRUM 

Plasma instabilities that lead to the formation of the collisionless 
shock are also responsible for the amplification of magnetic fields 
in the shock transition (Medvedev & Loeb 1999). Whether these 
small scale fields survive at large distances downstream from the 
shock is still an open question (Chang et al. 2008, Keshet et al. 
2009, Gruzinov 2008; Medvedev & Zakutnyaya 2009). This is, 
of course, very important, since the fields must remain strong on 
macroscopic scales for the synchrotron model to explain the after- 
glow observations (Rossi & Rees 2003, Waxman 2006). 

In this paper, we assume that the magnetic fields survive far 
downstream from the shock and facilitate the synchrotron and SSC 
radiation. In general, synchrotron dominates the X-ray and softer 
emission, while inverse Compton (IC) dominates in the 7-ray band 
(e.g., see Fan et al. 2008 and references therein). Here, we will 
neglect the IC component and only focus on the new features that 
appear in the synchrotron emission due to the thermal component 
in the electron distribution. 

We assume that a fraction eg ~ 10"- of the energy dissipated at 
the shock remains in the magnetic fields at macroscopic distances 
downstream from the shock. For a fluid moving with F » 1 just 
behind the shock (in the frame of the central engine), the energy 
density of the shocked fluid is e = 4T^ni,„mpC^ , where n,,,, is the 
number density of the external medium (e.g.. Sari et al. 1998). The 
comoving magnetic field strength is thus 



B- = SKege = 32;resF^n, 



(4) 



Using the standard expressions for synchrotron emission of 
ultrarelativistic electrons (Rybicki & Lightman 1979), we calculate 
the spectrum for various values of 6. In Fig. 2 we show a "slow 
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Figure 2. Synchrotron spectrum in the case of slow cooling for dilferent 
values of the fraction S of the energy in the power-law component of the 
electron distribution. The cooling frequency is taken to be = iO'^Vd, (vrf, 
is the characteristic synchrotron frequency). Self-absorption is not included 
in the calculation. 



, le-4 



le-5 



le-6 



le-7 



le- 



1 1 1 MMII| 


1 1 MIMI| 1 1 MIMI| 1 1 MIMI| 1 1 IIIMI| 1 1 Mllll| 1 


iiiiiii| 1 


1 Mllll| I 


Mill 














- - 5-1 \ 










— 5-1/3 \ 










— 5-0.1 










— 5-0.01 








1 


1 1 1 1 1 1 


1 1 


1 


^ 



e-4 le-3 0.01 0.1 



1 10 

v/v , 



100 le-l-3 le-H4 le-l-5 



Figure 3. Synchrotron Spectrum in the case of fast cooling for different 
values of the fraction of the energy in the nonthermal component S. The 
cooling frequency is taken to be = O.Olvrf,. The different particle distri- 
butions affect the spectral shape at the characteristic synchrotron frequency 
and the normalization of the high-energy power-law emission. 



cooling" case, where the cooling frequency lies above the charac- 
teristic synchrotron frequency. A reverse situation, corresponding 
to the "fast-cooling" case, is shown in Fig. 3. The cooling frequency 
is defined as the peak synchrotron frequency for electrons whose 
synchrotron cooling timescale equals the expansion timescale of 
the system. 

The spectra for different S are very similar below the charac- 
teristic frequency v^, = ^ylij^Vcyc, where v„,c = eB/lmn^c is the cy- 
clotron frequency. The Maxwellian declines fast enough for y <c 
that the emission at lower frequencies comes from particles close 
to the peak of the thermal distribution, which is close to the low 
energy cutoff of the pure nonthermal power law for (5 = 1 . Here we 
ignore the self absorption which appears in the low-frequency part 
of the spectrum (typically, radio). We estimate, however, that vary- 
ing 6 introduces a weak variation in the self-absorption frequency. 
Thus, we do not discuss the self-absorption further. 



The cooling frequency (or, equivalently, the Lorentz factor of 
electrons that emit at the cooling frequency) does not depend on 
the shape of the particle distribution. The difference in spectra due 
to particle distributions with different S is clearly seen around the 
characteristic frequency Vch ■ For 6=1 the characteristic frequency 
appears as a break in the spectrum. For 6 I, above the character- 
istic frequency (determined by the temperature of the thermal com- 
ponent) there is a sharp decline of the emission that is followed by 
hardening at higher energies. The sharp decline is not described by 
a power-law spectrum. The hardening corresponds to the emission 
coming from the nonthermal component of the electron distribu- 
tion. 

This particular shape around the characteristic frequency is 
reflected in the afterglow lightcurves. When the characteristic fre- 
quency crosses the observer band, new afterglow features appear as 
we discuss in the following section. 



4 AFTERGLOW LIGHTCURVES 

We consider a blastwave that decelerates due to the accumulation 
of mass from the circumburst medium. The bulk Lorentz factor of 
the blastwave depends on the total energy E, distance R and the 
external medium density (and profile) 



TiR) 



(5) 



where M„, = AnR'-rig^mpdR' is the total mass accumulated 
from the external medium at distance R. 

We focus on two types of external medium: a constant density 
medium, n,,,, =const, and a wind-like profile n„-, oc IjR-. Note, 
that we keep E constant, i.e., we ignore the energy injection that 
appears to be needed to explain the shallow decay segments of the 
X-ray lightcurves (Zhang et al. 2006; Nousek et al. 2006). Such 
inclusion is straightforward but does not affect our main points of 
discussion. 

Although the modeling of the afterglow emission involves sev- 
eral parameters, most of them have been extensively explored in 
the literature. Here we focus on the new effects on the lightcurves 
from varying the fraction 6 of the energy in the nonthermal com- 
ponent. We compute the lightcurves for E = 10^^ erg, = 0.3, 
€b = 0.01, and p = 2.5, which we will refer to as reference val- 
ues of the parameters. For constant density of the external medium 
we use riexi = 1 cm"^ as a reference value. For the wind case, the 
density profile is calculated assuming a spherical stelar wind of 
M = 10"^ Mo/year and velocity v,i. = 10* cm/sec expected from 
a Wolf-Rayet progenitor (see Li & Chevalier 1999). 

In computing the afterglow lightcurves, we use the standard 
approach described in Sari et al. (1998) and assume a burst at lumi- 
nosity distance d = 10"** cm. In Figs. 4 and 5 we show the resulting 
lightcurves in the optical and in Fig. 6 the lightcurve at 1 keV. 

The lightcurves are very similar for different d until the charac- 
teristic frequency crosses the observed band. For (5=1 the crossing 
appears as a single break in the lightcurve. On the other hand, for 
5 < 0. 1 there is a steep decline of the lightcurve followed by a break 
that leads to a more shallow decline. During the steep decline the 
spectrum shows a characteristic curvature (see next section). The 
time at which the steep decline is interrupted by a more shallow one 
occurs when the nonthermal electrons dominate the emission in the 
observed band. This shape of the lightcurve is a robust expectation 
from a mixed thermal-nonthermal electron distribution for (5 < 0.1. 
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Figure 4. Synchrotron lightcurve in tlie optical band for tlie reference values 
of parameters and constant density external medium. The various curves 
correspond to different values of S. The jet-spreading effects (resulting in a 
break in the afterglow lightcurve, so called jet break) are not included in the 
calculation. 
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Figure 6. X-ray lightcurve for the reference values of parameters and con- 
stant density external medium. A wind-like medium results in very similar 
lightcurves (not shown for clarity of the plot). The jet-spreading effects (re- 
sulting in the jet break) are not included in the calculation. 
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Figure 5. Same as Fig. 4 but for the wind-like external medium. 

4.1 Application to the early, steep X-ray decay 

Very steep decline in the early X-ray afterglow has been frequently 
observed hundreds of seconds after the GRB (Nousek et al. 2006). 
The steep decays were initially attributed to ofF-axis GRB emission 
(high-latitude model; Kumar & Panaitescu 2000). This interpreta- 
tion predicts, however, a specific relation between the temporal and 
spectral power-law indexes during the very steep decay that is not 
in agreement with that observed in the majority of the afterglows 
(O'Brien et al. 2006fl 

Zhang et al. (2007) did a systematic study of 44 steeply decay- 
ing X-ray afterglows for which time resolved spectra are available. 
They found that 1 1 X-ray tails did not show significant spectral evo- 
lution with time and are compatible with the high-latitude emission 

' Zhang et al. (2009) modified the high-latitude model by allowing for a 
non-power-law spectrum upon the cessation of the prompt emission phase. 
This model can account for the observed spectral evolution during the steep 
decline of the X-ray afterglows if at the end of the prompt emission the 
spectrum is steepening just above the XRT band. 



model. The rest show a clear spectral evolution not expected from 
the high-latitude emission model. Out of those 33 afterglows, 16 
have smooth tails (i.e., no flaring activity). The spectral index of 
the smooth tails showed a characteristic steepening with time in 
the X-rays from spectral index { fy ~ v^^) p ~ 0.5 to j8 ~ 2 (see 
Fig. 2 of Zhang et al. 2007). Moreover, after the very steep decay 
is over, the spectrum hardens top ~ 1.2. 

All these temporal and spectral features appear naturally when 
the characteristic frequency from a mixed distribution crosses the 
X-ray band for 6 ~ 0.1 and depend weakly on the parameters used 
to model the afterglow emission. 

Since Vrf, ~ ry^^^.^jVcc, the characteristic frequency as a func- 
tion of observer's time is (ignoring redshift corrections) 

vc, ~ 10' elo^e'J^El'^rJi' eV. (6) 

For deriving the last expression, we have used Eqs. (4), (5) and that 
the observer tim41?o*.i ~ R/4r^c. Note that Eq. (6) does not depend 
on the density of the external medium and approximately applies 
for both constant density and wind-like media. The characteristic 
frequency crosses the X-ray band marking the onset of the steep 
decay time t^d at 

t,,d ~ 300e%4%El'^ sec (7) 

after the burst (with the exact numerical value depending on S), 
leading to a steep decline of the lightcurve if (5 < 0. 1 (see Fig. 6). 

During the steep decline there is a characteristic softening of 
the spectrum, which hardens again after the the decay enters the 
power-law phase (see Fig. 7). This spectral evolution with time de- 
pends only on the fraction 6 and not on the external medium. The 
6p, £b and E parameters mainly affect when the steep decay (accom- 
panied by the spectral evolution) takes place but not the range over 
which the spectral index varies. 

In Fig. 8 we overplot our model predictions with afterglow 
observations. We select afterglows with good quality of data from 
the Zhang et al. (2007) sample that show the typical hard-soft-hard 
evolution. The observed steepening of the spectrum from p ~ 0.5 

^ The expression for t^bs is applicable within a factor of 2 for both wind 
and ISM external media. 
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Figure 7. Evolution of the spectral index /? (where fy ~ v~^) of the ~1 
keV emission as a function of time for a constant density external medium 
and for dilferent values of S. A wind profile leads to very similar results. 
The rest of the parameters are kept to their reference values. The spectral 
softening takes place simultaneously with the steep decay of the X-ray flux 
(see Fig. 6). 



(characteristic of fast cooling) to/? ~ 2 appears naturally for (5 ~ 0. 1 
while the hardening to /? ~ 1 after the steep decline is over cor- 
responds to the high-energy fast-cooling /?/2 segment of the syn- 
chrotron spectrurr|j. 

Energy injection in the blastwave may be needed to explain 
the shallow decay of the lightcurve observed after the steep decay 
phase (Zhang et al. 2006; Nousek et al. 2006). While the energy 
injection (that is not included in our calculation) affects the tem- 
poral evolution of the lightcurves, it does not affect the hard-soft- 
hard evolution of the spectrum that uniquely reflects properties of 
the particle distribution. We can, therefore, conclude that in our in- 
terpretation of the steep decline, 8 is constrained to be ~ 0.1 in 
agreement with recent PIC simulations. 



5 CONNECTIONS TO THE PROMPT GRB EMISSION 

So far we have applied a mixed distribution of accelerated elec- 
trons in relativistic shocks to the afterglow emission. figuThe dissi- 
pative process responsible for the prompt GRB emission is more 
uncertain, with internal shocks (Paczynski & Xu 1994; Rees & 
Meszaros 1994) and magnetic dissipation (Usov 1992; Thompson 
1994; Spruit et al. 2001; Lyutikov & Blandford 2003) being the 
leading candidates. The prompt emission may be the result of syn- 
chrotron (Katz 1994), photospheric emission ( Eichler & Levinson 
2000; Meszaros & Rees 2000; Pe'er, Meszaros & Rees 2006; Gian- 
nios 2006) or a combination of the photospheric (that dominates the 
~1 MeV band) and synchrotron, SSC emission (contributing from 
the optical to the ~ GeV band) because of gradual energy release 
(Giannios 2008). 

The internal shock model predicts mildly relativistic colli- 
sions. Furthermore, the colliding ejecta may be substantially mag- 
netized. Mildly relativistic collisions of (potentially) magnetized 
plasma have not been studied in the same detail by PIC simulations. 



^ Note that one of the models shown in Fig. 8 (see yellow line) shows an 
additional late-time raise of the spectral index from /? = 0.55 to /? = 1.1 
because of the cooling break crossing the X-ray band. 



Figure 8. The evolution of the spectral index fi in the X-rays as a function 
of time for a constant density external medium as predicted by our model. 
The dashed, red, black and blue lines are the same as in Fig. 7. The green 
line shows the model predictions for 5 = 0.1, eg = 10"^^, p = 2.2 and the 
yellow line for S = 0.1, E = 10^^ erg, eg = 10"^, and p = 2.2. The rest of 
the parameters are kept to their reference values. Overplotted is the spectral 
evolution of the GRBs 050814, 060428B and 060614 observed with XRT 
and analyzed in Zhang et al. (2007). 
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Figure 9. The fast-cooling synchi'otron spectrum in Ef(E) representation 
(arbitrary units) for a shell collision with parameters discussed in the text. 
A clear thermal-like bump appears when 5 < 0. 1 , followed by the high- 
energy emission with spectral slope of p/2. The low-frequency slope is 
that of fast-cooling electrons, i.e., with spectral slope fS = 1/2. 



Keeping all these caveats in mind, we explore the possibility that 
the prompt emission is produced due to synchrotron emission from 
internal shocks and that the distribution of the downstream parti- 
cles is a mixed, thermal-nonthermal one. The effect to the prompt 
emission coming from an injected strong Maxwellian component in 
the electron distribution has also been discussed in Baring & Braby 
(2004) (see also Pe'er et al. (2006)). 

As an example, we consider two shells ejected with a time dif- 
ference St = 0. l&_i sec. For the ratio of their bulk Lorentz factors 
r2/ri of a few, the shells collide at distance Rjs ~ T-^i^cSt/S, where 
F,;, is the Lorentz factor of shocked plasma (e.g., Bosnjak et al. 
2009). One can check that for eg = 0.01, = 0.3, F,,, = 100, the 
characteristic energy of the synchrotron emission is E^, ~ 40 keV 
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Fig. 9 shows the fast-cooUng synchrotron spectrum for differ- 
ent 6 for the aforementioned internal shock parameters. For the pure 
power-law electron distribution {6 = 1), the emission spectrum ap- 
pears as smoothly connected power laws. A bump at the peak of the 
Ef(E) spectrum is pronounced for (5<0. 1 in addition to the low- and 
high-energy power laws. A similar feature has been found by Bar- 
ing & Braby (2004) for a pronounced Maxwellian component of 
the downstream electrons (although they adopted a different shape 
for the particle distribution). Note also that the peak of the spectrum 
increases by a factor of several for 5 < 0. 1 . 

There is evidence for "thermal-like" excess components in a 
number of bursts (Ryde 2005; Ryde & Pe'er 2009f|. Furthermore, 
a fraction of the bursts show a very steep decline above the peak 
of the spectrum. These bursts might be connected with very weak 
nonthermal acceleration in the shocks. The latter is, for example, 
possible if the colliding ejecta are pre-magnetized with large scale 
toroidal fields (Sironi & Spitkovsky 2009). On the other hand, a 
sample of bright bursts does not show any evidence for a significant 
Maxwellian component in the electron distribution when their spec- 
trum is fitted with a synchrotron model (Baring & Braby 2004). The 
strength of the bump can be used to constrain 6 in the synchrotron- 
internal shock interpretation of the GRB emission. 

The spectral slope bellow the peak of the synchrotron spec- 
trum is that of fast cooling electrons, i.e. fj = 1/2. This slope is too 
soft in comparison to that observed in the majority of the bursts. 
This is a well-known problem of the synchrotron interpretation of 
the prompt GRB emission and is not "cured" by the presence of 
a strong Maxwellian component in the electron distribution. Syn- 
chrotron spectra with spectral slope j8 ~ are expected bellow 
the synchrotron peak provided that the inverse Compton cooling 
takes place in the Klein Nishina regime (Derishev, Kocharovsky & 
Kocharovsky 2001; Wang et al. 2009). Harder spectra can be re- 
sult of photospheric emission and/or heating of the electrons over 
a timescale much longer than the cooling timescale (slow cooling 
model; Ghisellini & Celotti 1999; Stern & Poutanen 2004; Pe'er et 
al. 2006; Giannios & Spruit 2007; Vurm & Poutanen 2009). 



6 DISCUSSION/CONCLUSIONS 

Recent PIC simulations of relativistic shocks show that the electron 
distribution forms a thermal component downstream of the shock 
that receives the largest fraction of the dissipated energy in addi- 
tion to the nonthermal, power-law component (Spitkovsky 2008a,b; 
Martins et al. 2009). Here, we considered the effect of the ther- 
mal component on the afterglow spectra and lightcurves and on the 
prompt GRB emission. 

A strong Maxwellian component introduces new phe- 
nomenology when the characteristic frequency of the synchrotron 
spectrum crosses the observed band. Instead of a break in the 
lightcurve predicted by the pure power-law model for the electron 
distribution (e.g.. Sari et al. 1998), we find a steep temporal decline 
followed by a break and a more shallow decay (see Figs. 4, 5, 6). 
During these phases the spectrum shows a characteristic hard-soft- 
hard evolution (Fig. 7). The steep decline appears in the X-rays at 
~hundreds of seconds after the burst independently of the external 
medium density and profile (see Eq. (7)). 

Note, however, that the Ryde (2005) analysis tits the spectra with a black- 
body and a single power-law component while synchrotron emission from 
the mixed electron distribution results in a bump surrounded by power-law 
high- and low-energy emission with different spectral slopes. 



The very steep decay observed in the early X-ray lightcurves 
has spectral and temporal properties well studied thanks to the XRT 
on board of Swift (Nousek et al. 2006; Zhang et al. 2007). In the 
majority of bursts, the steep decay shows a characteristic "hard- 
soft-hard" spectral evolution that is in agreement to that expected 
from our model when the thermal component in the electron dis- 
tribution contains ~10 times more energy than the nonthermal one 
(see Fig. 8; in agreement with recent PIC simulations). 

Similar steep decline and spectral evolution is expected when 
the characteristic frequency crosses other bands. From Figs. 4, 5 
one can see that the steep decay appears at around 1 day after the 
burst in the optical band (or slightly later if there is energy injec- 
tion). During the steep decline, the spectrum is predicted to devi- 
ate from a power-law. There are several optical afterglows showing 
breaks and steep decline at ~ 1 day after the burst. Unfortunately, 
the time coverage of the optical lightcurves is often sparse and there 
are other potential physical sources for breaks (end of the energy 
injection, jet break) that complicate the interpretation of the obser- 
vations (Liang et al. 2008). It remains to be seen if optical observa- 
tions support the mixed electron distribution that we propose. 

Although the self-Compton emission is not included in our 
calculations, we expect a similar steep decay signature when the 
Comptonised component of the characteristic synchrotron fre- 
quency i/Jj ~ yl,j^Vci, crosses the observed band. Depending on pa- 
rameters, the steep decay can take place at ~1000 sec in the ~GeV 
band that is now accessible to observations thanks to the FERMI 
mission (see, e.g.. Fan et al. 2008). 

Provided that the prompt GRB emission is the result of inter- 
nal shocks, similar considerations for the distribution of acceler- 
ated particles can be applied to the prompt GRB itself. If the large 
fraction of the energy goes into the Maxwellian component, the 
synchrotron peak of the prompt emission spectrum should exhibit 
a "bump" in addition to the typical fast-cooling synchrotron spec- 
trum from nonthermal (power-law) electrons that may be common 
in GRBs (Ryde 2005), but not universal (Baring & Braby 2004). 
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